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Optical beams with topological singularities have a Schmidt decomposition. Hence, they display features typ¬ 
ically associated with bipartite quantum systems; in particular, these classical beams can exhibit entanglement. 
This classical entanglement can be quantified by a Bell inequality formulated in terms of Wigner functions. We 
experimentally demonstrate the violation of this inequality for Laguerre-Gauss (LG) beams and confirm that 
the violation increases with increasing orbital angular momentum. Our measurements yield negativity of the 
Wigner function at the origin for LGio beams, whereas for LG20 we always get a positive value. 
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I. INTRODUCTION 


Entanglement is usually presented as one of the weirdest 
features of quantum theory that depart strongly from our com¬ 
mon sense • Since the seminal work of Einstein, Podolsky, 
and Rosen (EPR) H, countless discussions on this subject 
have popped up . 

A major step in the right direction is due to Bell |0], who 
formulated the EPR dilemma in terms of an inequality which 
naturally led to a falsifiable prediction. Actually, it is common 
to use an alternative formulation, derived by Clauser, Horne, 
Shimony and Holt (CHSH) H, which is better suited for re¬ 
alistic experiments. 

The main stream of research lill settled the main concepts 
of this topic in the realm of quantum physics. However, in re¬ 
cent years a general consensus has been reached on the fact 
that entanglement is not necessarily a signature of the quan¬ 
tumness of a system. Actually, as aptly remarked in RefE. JH], 
one should distinguish between two types of entanglement: 
between spatially separated systems (inter-system entangle¬ 
ment) and between different degrees of freedom of a single 
system (intra-system entanglement). Inter-system entangle¬ 
ment occurs only in truly quantum systems and may yield to 
nonlocal statistical correlations. Conversely, intra-system en¬ 
tanglement may also appear in classical systems and cannot 
generate nonlocal correlations JH] ; for this reason, it is often 
dubbed as “classical entanglement”. Since its introduction by 
Spreeuw 111^ . this notion has been employed in a variety of 
contexts ifllll . 

Classical entanglement has allowed to test Bell inequalities 
with classical wave fields. The physical significance of this vi¬ 
olation is not linked to quantum nonlocality, but rather points 
to the impossibility of constructing such a beam using other 
beams with uncoupled degrees of freedom. However, all the 
experiments conducted thus far to observe this violation have 
involved only discrete variables, such as spin and beam path 
of single neutron s 111^ . polarization and transverse modes of 
a laser beam |[13 | - [l^ . different transverse modes propagating 
in multimode waveguides 10, polarization of two classical 
fields with different frequencies 111^ . orbital angular momen¬ 
tum (HOjIziI], and polarization and spatial parity JH]. 


In this paper, we continue the analysis of this classical en¬ 
tanglement by focusing on the simple but engaging example 
of vortex beams. To this end, in Sec. m we revisit a decom¬ 
position of Laguerre-Gauss (LG) beams in the Hermite-Gauss 
(HG) basis that can be rightly interpreted as a Schmidt decom¬ 
position. This immediately suggests that many ideas ensuing 
from the quantum world may be applicable to these beams as 
well. In particular, in Sec. [Till we address the inseparability 
of the LG modes using a CHSH violation that we quantify 
in terms of the associated Wigner function. As this distribu¬ 
tion can be understood as a measure of the displaced parity, in 
Sec.[IV|we discuss an experimental realization which nicely 
agrees with the theoretical predictions. Einally, our conclu¬ 
sions are summarized in Sec. lYl 


II. OPTICAL VORTICES AND SCHMIDT 
DECOMPOSITION 


It is well known that the beam propagation along the z di¬ 
rection of a monocromatic scalar field of frequency co; i.e., 
E{rj) = (o (r)exp[—i{cot — kz)], is governed by the paraxial 
wave equation 
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with % — Xjln and A is the wavelength. Equation (12.1b is 
formally identical to the Schrodinger equation for a free parti¬ 
cle in two dimensions, with the obvious identifications t z, 
Xj/ 1 —y and h 1 —y A. 

Any optical beam can be thus expressed as a superposition 
of fundamental solutions of Eq. (12.1b . In Cartesian coordi¬ 
nates, a natural orthonormal set is given by the Hermite-Gauss 
(HG) modes: 
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where w is the beam waist, and are the Hermite polynomi¬ 
als. Note that we are restricting ourselves to the plane z = 0, 
since we are not interested here in the evolution. 

For cylindrical symmetry, it is convenient to use the set of 
Laguerre-Gauss (LG) modes, which contain optical vortices 
with topological singularities; they read 


function have been reported llsol - t^ . Note that W has the di¬ 
mensions of an intensity and it yields a description displaying 
both the position and the momentum (which in the paraxial 
approximation has the significance of a scaled angular coor¬ 
dinate) of the intensity of the wave field: in fact, one easily 
proves that 
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where l}p {x) are the generalized Laguerre polynomials. A 
word of caution seems to be in order: usually, these modes 
are presented in terms of two different indices: the azimuthal 
mode index i = m — n, which is a topological charge giving the 
number of 2;r-phase cycles around the mode circumference, 
and p = min(m,^) is the radial mode index, which is related 
to the number of radial nodes . However, the form (12.31) 
will be advantageous in what follows. 

The crucial observation is that the LG modes can be rep¬ 
resented as superpositions of HG modes, and viceversa. This 
can be compactly written down as in 

m+« 

LG^^(p,(p)= HG^^^_^j^(v,y) (2.4) 

k=0 

where the coefficients are 


t=0' 

(2.5) 

This looks exactly the same as a Schmidt decomposition for 
a bipartite quantum system. It is nothing but a particular way 
of expressing a vector in the tensor product of two inner prod¬ 
uct spaces 11^ . Alternatively, it can be seen as another form 
of the singular-value decomposition which identifies the 

maximal correlation directly. In quantum information, the 
Schmidt coefficients convey complete information of the 
entanglement na. Here, we intend to assess entanglement 
in LG beams via the violation of suitably formulated Bell in¬ 
equalities. 
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III. CHSH VIOLATION FOR LAGUERRE-GAUSS MODES 


j W (x,p) dp = /(x) = {E* (x)£(x), 

^ [w{x,p)dx = I{p) = {E*{p)E{p)), 
7t ^ 


(3.2) 


with 


E{,p) = -^ (E{x)&xp{ip-x/X)dx. (3.3) 
X J 

Thus, the marginals of the Wigner function are the intensity 
distributions in x or p space, respectively. 

The CHSH inequality can now be stated in terms of the 
Wigner function as 

B = ^\W{a,P) + W{a,l5') + W{a',l3)-W{a',l3')\<2, 

(3.4) 

where a = {x,px)/V2 and P = {y^py)/V2. This also follows 
from the work of Gisin iH, who formulated a Bell inequality 
for the set of observables with the property = t: as we 
shall see, the Wigner function appears as the average value of 
the parity, whose square is unity. Reference presents a 
detailed study of the violations of (13.4b . 

For the state LG^„, the normalized Wigner function can be 
written as 

Wk°{X,Px-,Y,Py) = ^ 4 exp(-46o) 

xL^[4(Go + e2)]i^«[4(Go-e2)], (3.5) 

where 

Q^=^{X^ + Y^ + Pl+Pi), Q2 =\{XPy-YPx), 

(3.6) 

and we have rescaled the variables as v i-^ (w/ V2)X and px ^ 
{pJlXjw) Px (and analogously for the y axis). Let us first look 
at the simple case of the mode LGio, which reduces to 


The traditional form of the CHSH inequality applies to di¬ 
chotomic discrete variables. For continuous variables, the sen¬ 
sible formulation is in terms of the Wigner function, which for 
a classical beam reads 

W(x,p) = ^ /c;Ve2'P-'/y£*(x-x')£(x+xO), (3.1) 
X J 

the angular brackets denoting statistical average. Although 
originally introduced to r^resent quantum mechanical phe¬ 
nomena in phase space 11^ . the Wigner distribution was es¬ 
tablished in optics to relate partial coherence with ra- 
diometry. Since then, a great number of applications of this 


wI^^{X,Px-,Y,Py) = T exp(-P| -P^-X^- Y^) 

x[{Px-Yf + iPY+Xf-l]. (3.7) 

The two measurement settings on one side are chosen to be 
a = (A = 0,Fx = 0) and a' = (A' =X,P^ = 0), and the cor¬ 
responding settings on the other side are j3 = (T = 0,/V = 0) 
and P' = (T' = 0,/^ = Pr) ll^ . for which the Bell sum is 

B = e~^y{P} - l)+e~^\x^ - 1) 

(3.8) 
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Upon maximization with respect to X and Py, we obtain 
the maximum Bell violation, |5max| — 2.17, which happens 
for the choices X ^ 0.45, Py 0.45 ll^ . For comparison, 
note that the maximum Bell violation in quantum mechanics 
through the Wigner function for the two-mode squeezed vac¬ 
uum state using similar settings is given by | 2.19 iH. 

The Bell violation may be further optimized by a more gen¬ 
eral choice of settings than those used here. For example, 
maximizing it with respect to the parameters a = (X^Px), 
a' = {X',Px), P = {Y,Py), P' = {Y'.Py), one obtains the abso¬ 
lute maximum Bell violation, |5max | = 2.24 and occurs for the 
choices X - -0.07, Px ^ 0.05, X' - 0.4, P^ - -0.26, Y - 
—0.05, Py —0.07, Y' 0.26, Py — 0.4. The violation also 
increases with higher orbital angular momentum. This in¬ 
crease with n is analogous to the enhancement of nonlocality 
in quantum mechanics for many-particle Greenberger-Horne- 
Zeilinger states (H]. 


IV. EXPERIMENTAL RESULTS 


We have carried a direct measurement of the Bell sums for 
optical beams with different amount of nonlocal correlations. 
To understand the measurement, we recall that the Wigner 
function in quantum optics is often regarded as the average 
of the displaced parity operator j^ . At the classical level, 
we can consider the field amplitudes (X, T) as vectors in the 
Hilbert space of complex-valued functions that are square in- 
tegrable over a transverse plane. In this space we define linear 
Hermitian operators 


X:^{XJ)^X^{XJ), 




(4.1) 

and analogous ones for the Y variable. Formally, these op¬ 
erators satisfy the canonical commutation relations [X^Px] = 
[Y,Py] = i. Therefore, the unitary parity operator is 


nxXnx = —X, tlxPxYlx =—Px , (4.2) 



FIG. 1. (Color online) Scheme of the Bell measurement. The abbre¬ 
viations are as follows: He-Ne: laser source, FC: fiber coupler, SMF: 
single mode fiber, CO: collimation optics, SLM: spatial light modu¬ 
lator, AS: aperture stop, BS: beam splitter, M1-M4: mirrors, CCD: 
camera 



FIG. 2. (Color online) Snapshots of the CCD camera for the state 
LGio at the four settings (X^Px'X^Py) indicated. The scans are nor¬ 
malized to the peak intensity among the measurements and the area 
of interest for the intensity integration is marked by a red circle. 


and changes ^(X,F) into ^(—X,F). The displacement oper¬ 
ators are 


D{X,Px) = txi^[i{PxX-XPx)\ . (4.3) 

Indeed, with these notations we have 
W{X,P,-Y,Py) = 

^{D{X,P,)nxD\X,P,)D{Y,Py)UYD^{Y,Py)). 

(4.4) 

Parity measurement can be, in turn, realized by a common- 
path interferometer with a Dove prism inserted into the opti¬ 
cal path El. In our setup, sketched in Fig.[T] the prism was 
substituted with an equivalent four-mirror Sagnac arrange¬ 
ment . The two copies of the input signal obtained after 
the input beam splitter are transformed by the mirrors so as 
to make one copy spatially inverted with respect to the other, 
prior to combining the beams together. The resulting interfer¬ 
ence pattern is detected by a CCD camera: Figure [2] shows 
snapshots of the camera for the state LGio at the four set¬ 
tings indicated. The total intensity witnessing parity of the 
measured beam is computed by spatial integration and this is 
proportional to the desired Wigner distribution sample after 
normalization to the overall intensity. 

The target signal beams were prepared with digital holo¬ 
grams created by a spatial light modulator (SLM), which mod¬ 
ulated a collimated output of a single mode fiber coupled to a 
He-Ne laser. We also included a 4/-system, with an aperture 
stop, to filter the unwanted diffraction orders produced by the 
SLM. To allow for a better fiexibility, all the necessary shifts 
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FIG. 4. (Color online) Experimental results for the beam LG 20 , 
presented in Fig. 3. The values of (X, 0; 0, Py) are also indicated. The 
Bell sum variation is significantly larger than in the case of LGiq. 
The plots on the right bottom panel are the amplitude and phase of 
LG20. 



a b c 


0.0, Px ^ 0.0, X' - -0.45, P^ - 0.0, Y - 0.0, Py - 
0.0, Y' 0.0, Py — —0.45 for the evaluation of the Bell 
sums. The theoretical values of the Bell sums for these are 
(-1.91,-2.15, —2.17), respectively. 





FIG. 3. (Color online) Experimental results for three different opti¬ 
cal beams: a) HGio, b) 0.4HGio + /O. 6 HG 01 , and c) LGiq. At the 
top, we plot {X,Px\Y,Py) at the values (X,Px\Y,Py) indicated 

for each one. The next plot shows the measured Bell sums, all re¬ 
ported with 75% and 25% quartile (orange boxes) and the minimal 
and maximal measured values (error bars). The theoretical values 
(—1.91,-2.15,-2.17) are the dots and the black bar is at |5| = 2, 
which delimites the classically entangled states. The theoretical am¬ 
plitude (top) and phase (bottom) distributions of the measured beams 
are plotted bellow the chart. 


Each measurement was repeated many times with slightly 
different readings, due to laser intensity instabilities and CCD 
noise. These effects manifest as measurement errors, which 
can be estimated from the sample statistics. As the parity 
measurement requires to normalize the total measured inten¬ 
sity of the interference pattern with respect to the input beam 
intensity, a separate reading of the input beam intensity was 
performed. For each optical beam, the mean value of the Bell 
sum is reported. The results are summarized in Fig. [3j The 
Bell correlations grow with the coupling between the basis 
HGio and HGoi modes, with statistically significant violation 
of CHSH inequality by the second and third beams, as theo¬ 
retically predicted. 

We also show the measured values of the Wigner function. 
For both, HGio and FGio modes, the values of (0,0; 0,0) 
are quite close to —1. For classical beams, ours is one of the 
few measurements on the negativity of the Wigner function, 
though it has to be anticipated from the corresponding results 
in quantum optics B. We note that very early, March and 
Wolf ll^ had constructed an example of a classical source 
which exhibited negative Wigner function. 


in the X, 7, Px, and Py variables were incorporated into the 
SFM, so that each Bell measurement was associated with a 
separate hologram. 

The measured beams were coherent superpositions of 
Hermite-Gaussian beams in the form a HGio + ibYLGoi with 
{a = l,b = 0}, {a = 0.4, b = 0.6} and {a = 0.5, b = 0.5}, 
respectively. The first and the third are thus a pure Hermite- 
Gaussian beam and a pure Faguerre-Gaussian vortex beam, 
respectively. For all the beams we used the settings X 


Finally, we have checked the violation of CHSH inequality 
for the beam FG 20 . A beam with higher topological charge 
is more sensitive to setup imperfections, hence the Bell sum 
variation is significantly larger than in the case of FGiq. On 
the other hand, as shown in Fig. [H the increasing of the Bell 
sum for higher orbital angular momentum is clearly demon¬ 
strated: the theoretical value for FG 20 is —2.24, which agrees 
pretty well with the experimental results, lli^] . Note that the 
Wigner function at the origin for the FG 20 beam is positive, as 
expected. 
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V. CONCLUDING REMARKS 


In short, we have presented an experimental study of nonlo¬ 
cal correlations in classical beams with topological singular¬ 
ities iH. These correlations between modes are manifested 
through the violation of a CHSH inequality, which we have 
detected via direct parity measurements. Such a violation is 
shown to increase with the value of orbital angular momentum 
of the beam. As a byproduct of our measurements, we obtain 
negativity of the Wigner function at certain points in phase 
space for the HGio and LGio beams. Note that this has im¬ 
plications for similar studies with electron beams, for which 
vortices have been reported 

Though entanglement here does not bear any paradoxical 
meaning, such as “spooky action on the distance”, it still rep¬ 
resents a potential resource for classical signal processing. It 
might be expected that future applications of quantum infor¬ 
mation processing can be tailored in terms of classical light: 


the research presented in this work explores one of those op¬ 
tions. 

Furthermore, our results are relevant not only for a correct 
understanding of “classical entanglement”, but also for bring¬ 
ing out different statistical features of the optical beams, since 
it provides an alternative paradigm to the well developed op¬ 
tical coherence theory. 
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